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ABSTRACT 

The observed mass function for all known extrasolar giant planets (EGPs) 
varies approximately as M~ l for mass M between ~ 0.2 Jupiter masses (Mj) 
and ~ 5 Mj. In order to study evaporation effects for highly-irradiated EGPs in 
this mass range, we have constructed an observational mass function for a subset 
of EGPs in the same mass range but with orbital radii < 0.07 AU. Surprisingly, 
the mass function for such highly-irradiated EGPs agrees quantitatively with 
the M~ l law, implying that the mass function for EGPs is preserved despite 
migration to small orbital radii. Unless there is a remarkable compensation of 
mass-dependent orbital migration for mass-dependent evaporation, this result 
places a constraint on orbital migration models and rules out the most extreme 
mass loss rates in the literature. A theory that predicts more moderate mass loss 
gives a mass function that is closer to observed statistics but still disagrees for 
M < 1 Mj. 



Subject headings: molecular processes — planetary systems 
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Introduction 



The rate of mass loss during evolution of highly-irradiated EGPs has important impli- 
cations for our understanding of the origin of these unanticipated objects. The principal 
purpose of this Letter is to show that enough such objects have now been detected to permit 
construction of a mass function for highl y-irradiated EGPs, i n the mass range where strong 
evaporation effects had been predicted (IBaraffe et al.l 120041 ) . The derived mass function 
shows no evidence for mass loss; its resemblance to the mass function for all observed EGPs 
may instead be a constraint on migration theories. 



Our analysis (IHubbard et al.ll2006l ) is based on a single dimensionless parameter e which 
represents the efficiency of the mass loss process: 



<t>E 



B,mod 



where $ is the flux of escaping molecules (number per unit area per unit time), E B ^ mod 
is the gravitational binding energy of a molecule to the planet, as modified by the tidal 
potential, and S 1 * is the bolometric stellar flux (energy per unit area per u nit time) at the 
planet 's orbital distance, for a main-sequence star of solar mass and age. IHubbard et al. 
( 120061 ) calculate the EGP mass-loss rate as a function of time by coupling this parameter to 
evolutionary models for EGPs orbiting solar-mass stars at orbital radii in the range 0.023 



(20 04J) c an be reproduced by setting 



AU to 0.057 AU. The calculations of IBaraffe et al. 
e = 10~ 4 . A model scaled to the calcula tions of IWatson et al.l (1198ll ). as well as to the lower 



limit considered b y IBaraffe et al. (12006 ). corresponds to e = 10~ 6 . In this paper, we use the 
suite of models of IHubbard et al.l (120061 ) to investigate the cumulative effect of mass loss on 
an initial mass function for highly-irradiated EGPs. 



2. Initial Mass Function 



We assume an initial mass fun ction f(M, 0) defined as follows. The database of EGPs 
known in 2005 (IMarcy et al.l 120051 ) has a mass function (the mass M is multiplied in most 
cases by the unknown sine of the orbit inclination i) corresponding to 



f(M, 0) = dN/dM oc M~ 



(2) 



Our independent fit to a slightly differe nt database , the t able of M sin % and orbital semimajor 



Schneider 



jgood), n mid- 2006, gives a s imilar result 



axis a for reported EGPs as given by 
dN/dM oc M" 119 . As discussed bv Burrows et al.l (1200 lh and iMarcv et al.l (120051 ). when 



one corrects the underlying function dN/dM for random orbit inclinations, the index is 
essentially unchanged, and for an index of —1 it is exactly unchanged. 
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In this paper, we define the initial mass function (IMF) for highly-irradiated EGPs 
somewhat differently than the IMF for stars, the latter referring to the stellar birth function. 
Since the majority (~ 75%) of detected EGPs have orbital radii > 0.07 AU, they cannot 
have suffered significant mass loss from atmospheric escape over their lifetime. We refer 
to the latter class as "field" EGPs. Our hypothesis is that highly-irradiated EGPs are not 
formed in situ, but are remnants of field EGPs that have migrated inward to small orbital 
radii during the ~ 10 6 to 10 7 years that the star's initial planet-forming nebula persists. Let 
the IMF denote the initial mass function for highly-irradiated EGPs at age t ~ 10 6 to 10 7 
yr in a mass range ~ 0.2 Mj < M <~ 5 Mj, when these EGPs start their atmospheric 
erosion. This IMF may well differ from the observed mass function for field EGPs, since 
the migration mechanis m could have a mass d ependence. One prediction of the IMF for 
highly-irradiated EGPs (IDel Popolo et al.l 120051 ) suggests a depleti on of planets with M > 4 
Mj but that l ower- mass EGPs migrate inward readily (see also iTrilling et al.l (119981 ) and 
Trilling et al.l (120021 ) ) . In this paper we make the provisional assumption that the IMF has 
the same index as the mass function for field EGPs. 

Figure 1 shows our assumed IMF f(M, 0) (heavy solid line), arbitrarily normalized to 
unity at M = 1 Mj. We assume that EGPs born at a ~ several to several xlO AU are 
deposited at smaller orbital radii a ~ fewxl0~ 2 AU within the first fewxlO 6 years of the 
parent star's lifetime. To map the IMF onto a time- dependent ensemble of er oding EGPs, we 



fix the distance from the star at one of the four standard distances studied in Hubbard et al. 



( 120061 ). ranging from 0.023 to 0.057 AU. We then randomly choose an exoplanet of initial 
mass M n from the I MF an d allow it to lose m ass as a function of time t using either the 
Watson et al.l (119811 ) or the iBaraffe et al.l (120041 ) prescription. 



It is necessary to impose a low-mass cutoff to the IMF. As shown by iHubbard et al. 



( 120061 ). hydrogen-rich EGPs with masses < 0.2 Mj and initial entropies corresponding to 
isolated EGPs at ~ 10 6 years of age, have such large radii that their atmospheres are tidally 
unbound (independent of the mass- loss rate), and so we restrict our analysis to Mq > 0.2 

Mj. 



3. Mass Loss Models 



We have synthe sized time-dependent m ass functions f(M, t) for EGPs at the four orbital 
radii investigated by lHubbard et al.l (120061 ). Figure 1 shows resulting mass functions at t — 5 
Gyr after simultaneous mass loss and evolution of the EGP. In all cases, the mass functions 
are normalized to unity at M — 1 Mj. For fixed total initial mass of the ensemble, all 
subsequent mass functions would plot below the IMF due to mass loss, with the strongest 
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deviations at the lowest mass. However, we renormalize the theoretical curves to 1 Mj for 
comparison with the observed mass function data points. 



3.1. Lammer Model 



We denote by "Lammer Mod el" the evaporation theory of lLammer et al.l (120031 ) as 
incorporated in the predictions of iBaraffe et al.l (12004 ) . We do not separately simulate 
the implications of more moderate mass -loss rates investigated in more recen t publications 
(IBaraffe et al.l 120061 : lAlibert et al.l 120061 ) because these are bounde d by the IBaraffe et al. 



( 2004 ) rates and the Watson model. In agreement with the findings of Hubbard et al. ( 2006h 



5-Gyr mass functions corresponding to the Lammer model show the strongest effect of mass 
loss. At the smallest orbital radius investigated, a = 0.023 AU, the predicted 5-Gyr mass 
function (solid curve) has a positive slope for the entire mass interval plotted in Fig. 1, 
while the IMF has a negative slope. This behavior occurs because mass loss biased to the 
lowest-mass EGPs rapidly depletes their larger initial numbers. Indeed, the mass function 
plotted is alm ost entirely populated by remnants of EGPs with initial masses Mo > 3 — 4 Mj 
(see Fig. 7 of iHubbard et al.l (120061 )). For the largest orbital radii investigated, a = 0.057 
AU, after 5 Gyr there is a peak in the mass distribution at M ~ 2.3 Mj. 



3.2. Watson Model 

In contrast, the Watson model predicts essentially no deviation of the mass function 
from the IMF for M > 1 Mj. However, as we see from the curves in Fig. 1, for M < 1 Mj 
the Watson model predicts significant deviations from the IMF, with a maximum at about 
one Saturn mass for a = 0.057 AU and at even higher masses for smaller orbital radii. 



Observed Mass Function 



The database from which we construct the mass function for highly-irradiated EGPs 
compri ses ~ 40 objects (Table 1). These objects were selected from the list of all reported 
EGPs ( ISchneider 1 120061 ) according to the criteria 0.2 Mj < Msini < 5 Mj and a < 0.07 
AU. We selected a central mass bin of width 0.3 Mj, according to 0.9 Mj < Msini < 1.2 
Mj, containing four objects (Table 1), with dN/dM to be normalized to unity at this bin. 
Bins were chosen to be of equal width to the central bin, except for the smallest-mass bin 
0.2 Mj < Msini < 0.3 Mj, which is of width 0.1 Mj, because of the cutoff at 0.2 Mj, and 
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for the largest-mass bins, 1.5 Mj < Msini < 3 Mj and 3.0 Mj < Msinz < 4.5 Mj, which 
are five times wider because of the small number of high-mass objects in the sample. 

The resulting mass function is plotted in Fig. 2, with error bars determined by Poisson 
statistics. Also shown in Fig. 2 is the average orbital radius (in AU) of objects in the bin 
(number to the left of the error bar), and the average age (in Gyr) of the host stars in the 
bin (number to the right of the error bar). 

It is remarkable th at the mass fu n ction for these highly-irradiated objects is in agreement 
with the field IMF (see 



Marcy et al.l (120001 ) for a discussion of selection bias). 



5. Mass Function With sini Factor 

For comparison of the mass-loss models with observed statistics, we require a trans- 
formation of the theoretical mass function to a function of the new independent variable 
Msinz. Let f(M,t) be the mass function of an eroded ensemble of EGPs after time t. We 
write 

f(M,t) = f(q/sni,t), (3) 

where the new independent variable q = Msini. Now holding the observed quantity q fixed, 
we average the ensemble of eroded EGPs over all solid angles of presentation of their orbits 
to the observer dfl, according to 

9(q,t) = I ^/(g/sinz,t) = j\nf(q/y/T=J?,t), (4) 

where dfi = sinidi. As would be expected, averaging over sin 2 moves the maxima of the 
mass functions to slightly smaller masses. The resulting functions g(q,t = 5 Gyr) for the 
Watson model (curves to the left of Msinz = 1 Mj) and the Lammer model (curves to the 
right of Msinz = 1 Mj) are plotted in Fig. 2. 



6. Conclusions 

Our results strongly suggest that the mechanism for transporting EGPs to small orbital 
radii does not significantly change their mass function in the mass range considered here, and 
that unless a remarkable compensation of evaporation for migration has occurred, we have no 
evidence for evaporation effects. It is difficult to envision such a compensation. For example, 
to fit the observed mass function using the Lammer evaporation model, we would need to 
start with an IMF with an index ~ —3.6, but the fit would be highly sensitive to the ages of 
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the objects and would rapidly change with time. Instead, a parsimonious interpretation of 
the data would be that mass loss does not play a dominant role in even the hottest EGPs 
with M > 0.2 Mj. If we take the observed numbers of EGPs in the two innermost mass bins 
to be statistically valid, even the more moderate Watson model still disagrees with the data, 
but certainly less decisively than does the Lammer model. However, our arguments seem to 
point in the direction of requiring the mass-loss parameter e to be even smaller than 10~ 6 . 
All published th eories give values higher than this number; for example iTian et al.l (120051 ) 
and lYellel (bood ) have 10~ 6 < e < 10" 4 . 



We must mention some caveats. Comparison of theory with the observational database is 
complicated by the fact that a majority of the observed objects in the bins 0.6 Mj < M sin i < 
1.2 Mj are transiting EGPs (indicated by an asterisk in Table 1), and thus the masses in 
these bins are true masses without the sini ambiguity. A more consistent mass function 
can be derived once we have enough transiting EGPs to construct a statistically significant 
mass function for them alone. Also, we have not in cluded "Neptunes" in our analysis, 



although some highly-ir r adiated "Neptune s" do exist (IButler et al.l 12004 ; iMcArthur et al. 



2004; 


Santos et al. 


2004; 


Lovis et al. 


2006) 



20061 ). Our assumed sharp mass cutoff at 0.2 Mj is 
actually an approximation to a somewhat fuzzy boundary whose details depend greatly on 
the exact model for tidal effects an d on the initial entro py of such low mass hydrogen- 
rich objects. Our previous analysis ( Hubbard et al. 20061 ) confirms that such "Nepti 



tunes 



cannot be Jupiter-like EGPs (i.e., predominantly hydrogen), but must have large fractions 
of elements with Z > 2. Whether nature could ever form hydrogen-rich low-mass (Neptune) 
EGPs is highly debatable. We eagerly await observations of transits by "Neptunes" . 

As more transiting EGPs are detected (and these will be predominantly highly-irradiated 
ones), a more definitive statistical test of orbital- migration and mass-loss theories will be 
possible. 



This study was supported in part by NASA Grant NAG5-13775 (PGG), NASA grant 
NNG04GL22G, and through the NASA Astrobiology Institute under Cooperative Agreement 
No. CAN-02-OSS-02 issued through the Office of Space Science. 
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Table 1. EGPs Used for Calculation of Mass Function 



Planet Name 


e 


(M sin i) /Mj 


a (AU) 


Star Age (Gyr) 


btar Mass (M© ) 






0.2 to < 0.3 Mj 






HD 76700 b 


0.13 


0.197 


0.049 


4.52 


1.0 


HD 88133 b 


0.11 


0.22 


0.047 


9.56 


1.2 


HD 168746 b 


0.081 


0.23 


0.065 


3.75 


0.92 


HD 46375 b 


0.063 


0.249 


0.0398 


4.96 


1.0 


HD 109749 b 


0.01 


0.28 


0.0635 


10.3 


1.2 






0.3 to < 0.6 Mj 






HD 149026 b* 





0.36 


0.042 


5.8 


1.3 


HD 63454 b 





0.385 


0.0363 


1.0 


0.8 


HD 83443 b 


0.012 


0.41 


0.0406 


2.94 


0.79 


HD 75289 b 


0.034 


0.42 


0.046 


4.96 


1.05 


HD 212301 b 





0.45 


0.0341 


5.9 


1.0 


51 Peg b 


0.013 


0.472 


0.0527 


6.6 


1.06 


HD 2638 b 


o 


0.477 


0.0436 


3.0 


0.93 


BD -10 3166 b 


0.07 


0.48 


0.046 


4.18 


1.1 


HD 102195 b 


0.06 


0.492 


0.0491 


2.4 


0.93 


HD 187123 b 


0.023 


0.528 


0.0426 


5.33 


1.06 


OGLE-TR-111 b* 





0.53 


0.047 






HAT-P-1 b* 




0.53 


0.055 


3.6 


1.12 


OGLE-TR-10 b* 





0.54 


0.04162 










0.6 to < 0.9 Mj 






v And b 


0.012 


0.69 


0.059 


2.41 


1.3 


HD 209458 b* 





0.69 


0.0474 


4.72 


1.01 


TrES-1 b* 





0.759 


0.0394 


2.41 


0.87 


HD 330075 b 





0.76 


0.043 


6.21 


0.95 


WASP-2 b* 




0.88 


0.0307 






WASP-1 b* 




0.89 


0.0382 
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Table 1 — Continued 



Planet Name e (M sin i)/Mj a (AU) Star Age (Gyr) Star Mass (M Q ) 



0.9 to < 1.2 Mj 



XO-1 b* ••• 0.9 0.0488 

HD 179949 b 0.022 0.98 0.0443 2.05 1.24 

HD 188753Ab 1.14 0.0446 

HD 189733 b* 1.15 0.0312 6.1 0.9 

OGLE-TR-113 b* 1.19 0.0306 ••• 1.35 



1.2 to < 1.5 Mj 



TrES-2 b* 1.28 0.0367 ••• 1.08 

OGLE-TR-56 b* 1.29 0.0225 2 1.17 

OGLE-TR-132 b* 1.32 0.0229 > 0.7 0.78 

HD 149143 b 1.33 0.0531 7.6 1.21 

HD 86081 b 0.006 1.49 0.0346 ••• 1.0 



1.5 to < 3 Mj 



HD 68988 b 0.1249 1.86 0.0704 6.78 1.2 

HD 73256 b 0.029 1.87 0.0371 0.83 1.05 

HD 118203 b 0.309 2.14 0.0703 4.6 1.23 



3.0 to < 4.5 Mj 



HIP 14810 b 0.148 3.84 0.0692 • • • 0.99 

t Boo b 0.023 3.9 0.046 2.52 1.3 



* transiting object 
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Fig. 1. — Predicted mass functions f(M, t) for EGPs suf f ering mass loss over 5 Gyr, evaluated 
for four different orbital radii, using the IWatson et al.l (Il98ll ) escape rate (left-hand curves 
with maxima at masses be tween 0.4 Mj and 0.9 Mj), and the two-orders-of-magnitude 
larger lLammer et al.l (120031 ) escape rate (right-hand curves with maxima at masses > 2.3 
Mj). The heavy curve (IMF) is the assumed mass distribution at the start of mass loss. All 
distributions are normalized to unity at M = lMj. See text for a more detailed discussion 
of this plot. 
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Fig. 2. — Predicted mass functions g(M smi,t) for EGPs suffering mass loss over 5 Gyr, 
evaluated for the four different orbital radii and the two mass-loss models, compared with 
data for highly-irradiated EGPs. All distributions (theory and data) are normalized to unity 
at Msini = 1 Mj. See text for a more detailed discussion of this plot. 



